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Bose-Einstein condensate (BEC) is considered under conditions of Feshbach resonance in two- 
atom collisions due to a coupling of atomic pair and resonant molecular states. The association of 
condensate atoms can form a molecular BEC, and the molecules can dissociate to pairs of entangled 
atoms in two-mode squeezed states. Both entanglement and squeezing can be applied to quantum 
measurements and information processing. 

The processes in the atom- molecule quantum gas are analyzed using two theoretical approaches. 
The mean-field one takes into account deactivating collisions of resonant molecules with other atoms 
and molecules, neglecting quantum fiuctuations. This method allows analysis of inhomogeneous 
systems, such as expanding BEC. The non-mean-field approach — the parametric approximation — 
takes into account both deactivation and quantum fluctuations. This method allows determination 
of optimal conditions for formation of molecular BEC and describes Bose-enhanced dissociation of 
molecular BEC, as well as entanglement and squeezing of the non-condensate atoms. 



Introduction 



Recent developments in the field of Bose-Einstein condensation of dilute atomic gases are described in review 
articles [H 0, [3, 13, 0] and a book Q . One of the last accomplishments in this field is the formation of molecular 
Bose-Einstein condensate (B EC) from an atomic one 0, H, H, UOj UH [3 ^'^'i from a quantum-degenerate Fermi gases 
[H, Hi [S [ilj [13) [H, 0j HJ] ■ These experiments used the effect of Feshbach resonance (see Ref. [l!]), appearing 
in multichannel scattering when the collision energy of an atomic pair lies in the vicinity of the energy of a bound 
(molecular) state in a closed channel. An analog of this effect — the confinement induced resonance - appears in 
scattering under tight cylindrical confinement (see Refs. [23l. [2^ [25j ) . A cooperative effect of the confinement and 
closed channels has been considered in Ref. [1^ . 

Feshbach resonance allows control of BEC properties by tuning the elastic scattering length, as has been proposed 
in Refs. [27l . [?8. . 291. The open and closed channels can be coupled by resonant optical fields [H, [2^ or by hyperfine 
interaction [27| , using the Zeeman effect for the energy detuning control. A coherent formation of molecules due to 
such coupling schemes has been considered in Refs. |22l [sol. Isil. [33. Issj . The experimental realization of the scattering 
length control ^34, 35, 36] has demonstrated a drastic condensate loss when the resonance was approached. This loss 
has been attributed to two loss mechanisms: a deactivation of the temporarily formed molecules in inelastic collisions 
with atoms and molecules (see Refs. 



22, 
4C, 
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37, Ml) and a dissociation of the molecules resulting in the formation of non- 



condensate atoms (see Refs. 



14311'). Secondary collisions of the relatively hot deactivation products 
can lead to an additional loss (see Refs. |37l. |44|. lisjl. 

Many-body effects on the formation of molecules can be described in a mean-field approach by a set of coupled 
Gross-Pitaevskii equations [11] for the atomic and molecular mean fields (see Ref. [121) ■ The coupled fields form 
a hybrid atom-molecule condensate. Deactivation losses can be introduced into this approach by adding imaginary 
non-linear terms. A derivation of such equations is p resented in Sec. |TT] below. The mean field approach has been 
used for description of atom- molecule oscillations [U [H, H^] , allowing analytical solutions for some time-dependent 
models [isl . l49j . This approach describes also the effects of self-trapping and soliton formation [s^, [5l| , vorticity [s^l , 
and expansion (ssj of hybrid atom- molecule condensates. 

The second loss mechanism — the dissociation of molecules onto non-condensate atoms — has been described by 
a two-body theory in Ref. [39| and introduced as additional loss terms into the many-body mean-field equations in 
Ref. [i^. Both these approaches forbid the dissociation in a backward sweep, when the molecular state crosses the 
atomic ones in a downward direction. This manner, proposed in Ref. jsoj , has been used in experiments 
However, the dissociation can proceed even in a backward sweep, and correct description of this process requires more 
comprehensive theories. On a two-body level it can be described as a so-called "counterintuitive" transition [53 |. 
Apart from this, the dissociation of molecular BEC demonstrates quantum many-body effects beyond the mean-field 
approach, such as Bose-enhancement [H,!!!]. The effects of deactivation and dissociation losses are non-additive (see 
Ref. [H). 

The atomic pairs produced by dissociation are formed in two-mode squeezed states which are entangled [ssf. 
Squeezed states are characterized by noise reduction, and can be applied in communications and measurements. Two- 
mode squeezed states have the property of relative number squeezing (see Ref. [57]), which is important for atomic 
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interferometry. Entangled states of a decomposable system cannot be expressed as a product of the component states, 
and can be used in quantum computing and communications (see Refs. (58l [59|.) An entanglement measure has been 
considered in Refs. [13, EH. Some aspects of the entanglement have been considered within two-body theory of 
dissociation of individual molecules in Ref. fE^ , but the analysis of many-body effects requires a theory treating the 
non-condensate atoms as quantum fluctuations. 

Several theoretical methods are available for a more correct treatment of the many-body non-mean-field effects. A 
numerical solution of stochastic differential equations has been used in this context in Refs. [4^.[63l.[6^ 1 165", '66', '67', 'el, 
[69j . The role of quantum statistics of the initial state has been investigated by this method in Refs7|70., .71. 72, 73|. 
A direct solution of the many-body Schrodinger [tII and Liouville von-Neumann equations [t^ has been applied only 
to a small number of atomic modes (maximum two, in Ref. fl^)- A special time- independent case of a single atomic 
mode allows an exact solution in terms of the algebraic Bethe ansatz [t^]- Apart from this, non-mean field effects of 
a periodic potential [t^ and of the directional emission of correlated atomic pairs 79] have been analyzed. 

The prominent many-body quantum effects, such as entanglement, squeezing, and Bose-enhancement, can de de- 
scribed by methods taking into account second-order correlations. In the Hartree-Fock-Bogoliubov formalism of 
Refs. j43l [80l |. the mean field equations are complemented by equations for the normal and anomalous densities of 
atomic fluctuations. The macroscopic quantum dynamics approach [81] considers the molecules as two-atom bound 
states, incorporating actual interatomic potentials. The system dynamics is described by an inflnite set of equations 
for non-commutative cummulants. The flrst-order cummulant approach, applied to atom-molecule systems in Refs. 
[82, 83, 84, 85, 86], keeps the cummulants corresponding to mean flelds and anomalous densities. This truncation 
allows to express the contribution of quantum fluctuations in terms of two-body transition matrix. The resulting 
dynamic equations are simpler then the Hartree-Fock-Bogoliubov ones, allowing analyses of spatially inhomogeneous 
systems. The truncation neglecting the cummulants corresponding to normal density is justified under conditions of 
majority of current experiments. However, it leads to the loss of some quantum effects, such as Bose-enhancement 
(see Sec. IHI G[) . These effects could be taken into account in higher-order cummulant approaches. Such approach 
incorporating third-order correlations has been developed in Ref. [l20l ] for pure atomic gas. 

The theoretical approaches of Refs. [i^, IH, [6^-[8^ did not take into account the effects of deactivating collisions. 
The present chapter describes the parametric approximation — a non-mcan-field formalism in which the damping 
effects are incorporated (see Sec. IHip . This approach has been applied to the case of a single atomic mode in 
Refs. [75!, IstI], where analytical solutions have been obtained (see also Sec. IHI Fl and Ref. [1^), to the formation of 
instabilities in pure atomic BEG in Ref. 89], and to multimode atom-molecule quantum gases in Refs. [53 . 55, 56]. 
When the deactivating collisions are neglected, the parametric approximation becomes equivalent to the Hartree-Fock- 
Bogoliubov approach of Refs. [1^, [13] (see Sec. IIIIEl where Hartree-Fock-Bogoliubov equations involving deactivation 
are also derived). 

This chapter presents applications of the mean field and parametric approximations to the description of condensate 
losses (Sec. lIVI) . formation of molecules (Sec. El, and production of entangled atoms (Sec. lVI|) . The results demonstrate 
the importance of both non-mean-field and dumping effects for a correct description of processes in atom-molecule 
quantum Bose gases. A different situation takes place in the case of molecular BEG formed from Fermi atoms, where 
the deactivation can be neglected due to Pauli blocking (see Refs. [13, [Ui [13 ) • For a theory of molecular formation 
from quantum-degenerate Fermi gases, which is not presented here, see Refs. [9^ [13, [IE]- Production of entangled 
fermionic pairs has been analyzed in Ref. [96| . 

A system of units in which Planck's constant is fi. = 1 is used below. 



I. SECOND-QUANTIZED DESCRIPTION OF HYBRID ATOM-MOLECULE GASES 

The effect of Feshbach resonance appears when the collision energy of a pair of similar atoms of type A in an open 
channel is close to the energy of a bound state A2 (m) in the closed channel (see Fig. [T|). The temporary formation 
and dissociation of the resonant (Feshbach) molecular state A2 (m) can be described by the reversible reaction 

A + A^A2(m). (I.l) 

Although a resonance can be related to several closely-spaced bound states, this chapter considers a case of a non- 
degenerate resonant state, well separated from other bound states. This state is described by the field annihilation 
operator tpm (r), while the atoms are described by the operator 4'a (r)- The Hamiltonian for the free fields can be 
written in the form 

i^frco = / d^r^l (r) Ha^a (r) + / d^Tyn^i, (r™) Hyni'^ (r™) , (1.2) 
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FIG. 1: Schematic description of channel potentials and molecular states. The dashed line marks the open channel threshold. 



where 

Ha = TT^P' + Va (r) + e„, H,n - -^P™ + (r^) (1.3) 
Zm Am 

arc the Hamiltonians for the noninteracting atoms and molecules in the representation of first quantization. Here m 

is the atom mass, Va (r) and Vm (jm) are, respectively, the atomic and molecular trap potentials as functions of the 
position r of the atomic (or r„i of the molecular) center of mass. The time-dependent Zeeman shift 



(1.4) 



of the atom in an external magnetic field B (t) is counted from half the energy of the resonant molecular state (which 
is fixed as the zero energy point). Here 



(1.5) 



fia and /im are the magnetic moments of the atom and the resonant molecule, respectively, and Bq is the resonance 
value of B in the free space (when Va = Vm = 0). 

Since the atoms and molecules are treated here as independent particles, the interaction responsible for the atom- 
molecule coupling [reaction (jLl[) ] can be written in the general form 



14= / d^rd^r'y^(r-r')^t 



tpa (r) tpa (r') , 



(1.6) 



in which the molecule preserves the position of the center of mass. Although the interaction is localized within a range 
of atomic size, negligibly small compared to the condensate size and the relevant de Broglie wavelengths, the use of 
a zero-range interaction would lead to a divergence in the ensuing calculations. Therefore we keep this interaction as 
a finite-range function Vh (r — r') of the interatomic distance. 

The resonance bound state A2 (m) is generally an excited rovibrational state, belonging to non-ground states of the 
fine and hyperfine structures (see Fig. [T]) . This state can be deactivated by an exoergic collision with a third atom of 
the condensate [H [H IH, Si , 



A2(m) + A2(d) + A, 



(L7) 



bringing the molecule down to a lower state A2 (d) (see Fig. [T]) , while releasing kinetic energy to the relative motion 
of the reaction products. Although the collision occurs with a vanishingly small kinetic energy, rates of such inelastic 
processes remain finite at near-zero energies [12, 97, 98]. A variant of this process, involving deactivation by a collision 
with another molecule (rather than an atom), of the type 



A2(m) + A2(m) A2 (d) + A2 (u) . 



(I- 
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would require a significant molecular density to be effective. The two molecules emerge in two states A2 (d) and A2 (u) 
(see Fig.[T]), where A2 (u) can be a bound molecular state above A2 (d), or a continuum state of a dissociating molecule. 
The reaction can take place as long as the corresponding internal energies Ed and obey the inequality Ed + Eu < 0. 
A particularly effective reaction of type (jl.8[) would occur in the near- resonant case, in which < Eu < \Ed\. A typical 
example, common in VV-relaxation, is that oi v + v ^ {v + 1) + {v — 1) (where v is the vibrational quantum number 
of the state m). In this example, the kinetic energy is provided by the vibrational anharmonicity. 

The products of reactions (|L7} and (|I.8|) are described by the annihilation operators i^d (r) and ipu (r) , respectively, 
for each of the states A2 (d) and A2 (u) . The Hamiltonians of free fields can be represented as 

Hd,u = J d^rrai^X^ (''™) (i^P™ + ^'''") '^'■"^ ' 

where the trap potentials and Zeeman shifts can be neglected compared to the transition energies Ed^u- The deacti- 
vating collisions p.7p and (|I.8|) are described, respectively, by the following interactions 

Vd= I d\d\„Jd (|r - r™|) V-l (r) V-^ (r„0 ( l±pl] v3„, ( l±pA (HQ) 



3 / ' V 3 



Vud = J dVid3r2dw(|ri -r2|)^i(ri)^t (^^4^) ^'^^ (^4^) ' ^^'^^^ 

where dd and dud are the respective interaction functions in first quantization. The elastic collisions between atoms 
and molecules are described by the interaction term 

^+ (r) (r) V-a (r) i^a (r) + %i ^ (^) ^t' (^) (^) 



Vei = J d^'r 

+t/am^^+ (r) ^+ (r)V'. (r) (r) 



a,a' — m,u,d (1-12) 



including terms proportional to the zero-momentum atom-atom, molecule- molecule, and atom-molecule elastic scat- 
tering lengths (cq, a„i, and Oam, respectively), 

47r 27r 37r 

mm m 

The different numerical factors in the numerators reflect the different reduced masses. 

Finally, the total Hamiltonian of the atom-molecule gas is represented as the sum of free field Hamiltonians and 
interactions related to all collision processes taken into account, 

H = Hi,,, + Vu + Vl +Y,[Hd + Vd + + {^-d + V!d) + Vei- (1.14) 

d u u,d 

The Hamiltonian was described above in the Schrodinger picture using the coordinate representation. The momen- 
tum representation is more suitable for homogeneous systems. It can also be used for an approximate analysis of 
inhomogeneous systems in the local density approximation. The operator equations of motion in the Heisenberg 
picture will be used below for the analysis of non-mean-field effects. The Heisenberg field operator in the momentum 
representation (Pj t) is related to the corresponding Schrodinger operator in the coordinate representation ipa (r) 
as 

(p, t) = (2^)-3/2 J d're-'P'^Cl-^^ (r) Cl, (1.15) 

where the time evolution operator Q satisfies the equation 

d ^ 

i—n = Hn. (1.16) 
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As a result, the transformed Hamiltonian can be expressed in terms of the field operators in the momentum repre- 
sentation as 



p 

2m 



^-1 (p,i)*,(p,t) 

^-1 (p,0*c {P..t) \ + 



Am 



En 



-Vk + v,l + Y.{h + vj) + E {Vu, + v^) , 



(1.17) 

(1.18) 
(1.19) 
(1.20) 

(1.21) 
(1.22) 
(1.23) 

The transformation of the elastic collision terms is not presented here as the latter can be neglected in the analysis 
of non-mean-field effects below, carried out in the momentum representation. 



where the interaction operators are represented as 

V„ = {2tt)-'/' J dVV4(p-p')*L(P + P',i)*a(p,t)*a(p',i) 
Vd = (2^)-'/' / d^'d'pdd'pmdd (Pd - 2p') *t (p'^ t) ¥^ (pd, t) 

X i>rn (Pm, t) *a (p' + Pd - Pm, t) 

Vud = (27r)"^/^ / d^p„,d^pdd^pudud (p« - Pd) ^1 (p«, <) (Pd' *) 

X *„i (pm, t) (p„ -f Pd - Pm, t) , 

with 

d/i {p) = y rf^pl4 (p) exp (-ipp/2) 
dd (p) = (27r)"'/' y d^dd (p) exp (-ipp/3) 
d«d (p) = (27r)~^/^ / d^pdud (p) exp (-ipp/2) . 



II. MEAN FIELD APPROXIMATION 



Trial function 



Mean field equations can be derived by a variational method (see Refs. [99|, |100[) using a trial function in the form 
of a coherent state. This approach is applicable to the description of macroscopically occupied states, such as the 
atomic and resonant molecular states in the system considered here. However, "dump" states — "hot" products of 
the deactivating reactions p.7[) and p.8[) — are formed in a wide energy spectrum, escape the trap very fast, and 
therefore have low occupation. The form of Eqs. (jl.lOp and p. lip demonstrates that the resulting atoms and molecules 
are formed in these reactions as correlated pairs. These reasons suggest the following choice of a trial function (see 
Ref. ;40]), 



1$) = 



1 



d^rd^r,n ^ '^d (r, r,n,t) -0+ (r) ipj (r^) 



d^ridVa ^ V'w (ri, ra, i) -0+ (ri) ipj (rj) \ipo, fm) 

u,d 



It contains a coherent state 



\ipo, ifm) = exp| / d^r[ipo (r, t) 0+ (r) -f- ip^ (r, t) ^+ (r)] }|0), 



(III) 



(II.2) 



formed by a product of exponential operators, involving atomic (^o) ^-nd molecular {(pm) condensate states, and 
operating on the vacuum state |0). The linear factor preceding it in Eq. pi.l[) includes the fields (pd {^,^2,t) and 
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(pud (ri, r2, t), which are the correlated dump states of the products in reactions (jl.7|) and (jl.8|) . respectively. Another 
constraint, that follows from the large energy difference between the dump states and the resonant state, is the 
condition 



d^r^*oir,t)ifa{r,rra,t) = (II.3) 

[see discussion following Eq. pi.lip for justification]. 
The trial wavefunction pi.ip can now be substituted into a variational functional (see Refs. [o^ Il00l |) 



($1$) 

In the mean-field calculations one can use the approximation of zero-range interaction Vh (r — r') = gS (r — r'), and 
represent the atom-molecule coupling Vh [see Eq. (II. 6[) ] in the simpler form 

Vh^gJ (r) (r) (r) . (II.5) 

The parameter g is related to integral of Vh (p) and to the peak value of dh (p) [see Eq. (|I.21|1 ] as 

g = [ d^pVh (p) = dh (0) . (II.6) 



The interactions responsible for the deactivating collisions (|I.7|) and (jl.8|) are kept as finite-range functions of the 
distance between the reaction products [see Eqs. (jl.lOp and p. lip ]. 

Neglecting terms of th e order of (/s^ and (/s^, and taking into account Eq. (jTOj, the use of a standard variational 
procedure (see Ref. |lOOj |) then leads to a set of coupled equations for the atomic {ipo) and molecular {(pm) condensate 
fields (or "wavefunctions" ) , as well as for the dump states {tpd and (pud)- 

ipa (r, t) = {Ha + Ua\^Po (r, t) [^ + Uam\Vm (r, t) (^0 (r, t) 

+25 VS (r, t) ipm (r, t) + Q (r, t) p*^ (r, t) (11.7a) 

iiprn (r,t) = [Hyyi + Uam\Va -t) [^ + Um\Vni (r, Vni (r, 

+gipl (r, t) + Q (r, t) p* (r, t) + Q,„ (r, t) (r, t) (II.7b) 

itpdiri,r2,t) = [^pl + ^pl-Ed] pd{ri,r2,t) 
\ Zm Am J 

+d.(|ri - r^D^o ( ^-^^A ^™ ( "^^A (n-7c) 



iCpud (ri,r2,t) = (^^Pi + - -E'" - <^"d (ri,r2,i) 

+d„rf(|ri-r2|)^^('^:i±^,A (II.7d) 



where 



ri + 2r2 



Q(r,<) = /'dVidV2^dS(|ri-r2|)^d(ri,r2,i)(5('r- 
•' d ^ 

(r, = 2 1 ST^d\2 J2 did (|ri - r2|) (r, i) V'.d (ri, r2, i) x 5 (^r - ^^4^) 



(IE 



(II.9) 



B. Dump state elimination 

The procedure used in Ref. [13] to eliminate the dump states is similar to the Weisskopf-Wigner method in the 
theory of spontaneous emission (see Ref. [10i|). Equation (|lL7c)) is of the form of a Schrodinger equation for two free 
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particles with a source (the last term in the right-hand side) . Such an equation can be solved by applying the Green's 
function method for free particles, with the result 



ifd (ri,r2,t) - 



(27r) 



dt' / (fPd?p(fp(fRdd [p) (R, t') ipo (R, t') e-'PP-'^^ 



X exp 



' 6m 4m / 3 



(11.10) 



Here R is the center-of-mass position of the three-atom system, p is the radius vector of the reaction products, and P, 
p are the corresponding momenta. Since (pm (R, t) and ipo (R, t) are condensate wavefunctions, the Fourier transform 
of their product [the integral over R in Eq. (|II.10[) ] vanishes if P > 1/6, where b ~ (mwtrap) is a characteristic size 
of the condensate and Wtrap is the trap frequency. Therefore P^/ (6m) < Wtrap is negligible compared to Ed ^ t^^trap- 
This fact allows us also to neglect the time dependence of tpm (R, t) and (po (R, t) in the integration over t', and thus 
obtain the simplified expression 



ifd (ri,r2,t) 



1 



(27r)3 



ri + 2v2 



t (y9„ 



n + 2v2 



3 ,3 dd (p) exp [ip (ra - ri - p)] 



3p^ / (im) - Ed - iO 



The atom-molecule pair is thus formed with a momentum of relative motion 



Pd = ^imEd/3. 



(11.11) 



(11.12) 



The function ipd (ri, r2,t) is a rapidly oscillating function of the coordinates, and therefore condition (jll.3|l is justified. 
Substituting Eqs. (|ILTT|) and (|L22|) into Eqs. flLS)) one obtains 



where 



Q (r, i) = - (Sa + ika/2) ipo (r, t) ip„-, (r, t) , 

\dd [ip) ? 
3p^/ (Am) - Ed-iO' 



(11.13) 



(11.14) 



Using the well-known identity {x — iQ) ^ = Vx ^ -^ittS {x), where V denotes the Cauchy principal part of the integral, 
allows us to obtain explicit expressions for ka and 5a, 



in^m \ - I ~ , 



Sa 



IQirm 



Y^v [dp 



p'\dd (3p)| 
P'-Pl 



A similar analysis, starting from Eq. (jll.7dp . gives 

Qm (r, t) = ~ {5m + ikm) ipl^ (r, t) 

where 

kjn = STT^m ^ Pud dud (2P„d) 



167rm V [ dp^^ 
J P 



p'\dud{2p)\ 



2 „ p2 



(11.15) 
(11.16) 

(11.17) 

(11.18) 
(11.19) 
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Pud = V2m (Eu + Ed), and dud is defined by Eq. p.23|) . Substituting Eqs. (|II.13I) and (III.17|) into Eqs. (|II.7ap and 
(|II.7b[l one finally obtains a pair of coupled Gross-Pitaevskii equations (see Ref. 46]) 

iipQ = (^Ha + Ua\ipo\'^ + C^amlVmP) ^0 + 2g*ip^ipm - {Sa + ika/2) \iPm\'^ipo (II.20a) 

- [{5a + ika/2) |(/3oP + {6,n + ikm) Vm- (II.20b) 



The parameters 5a, ka, 5„i, and km, which are expressed in terms of dd and dud [see Eqs. pi.isp . pi.l6p . pi.lSp . 
and pi.l9|) ]. describe the shift and the width of the resonance due to the deactivating collisions with atoms and 
molecules, respectively. The parameters ka and km are the corresponding rate constants. The shifts 5a and Sm can 
be incorporated in the interactions Uam and Um, respectively. 

In the case of a time-independent magnetic field and large resonant detuning, and neglecting the decay described by 
the imaginary terms, Eqs. pi.20a|l and pl.20bp can be reduced to a single Gross-Pitaevskii equation with an effective 
scattering length [1 — A/ (_B — Bq)]. The phenomenological resonance strength A can be measured in experiments 
or obtained from two-body scattering calculations. It is related to the atom- molecule coupling constant g of Eq. pi.5[) 
as 

= 2TT\aa\fiA/m, (1121) 

from which the value of g can be deduced, given A. 

The atom-molecule coupling terms in Eqs. pi.20a[) and pl.20bp are proportional to the atomic mean field ipo. 
Therefore the dissociation of molecules into the atomic vacuum [ipo — 0) is forbidden. In the absence of molecule- 
molecule collisions {k„i = 0) a pure molecular BEG is a stationary solution of the mean-field equations. This solution 
is, nevertheless, unstable, and a proper description of the dissociation requires the use of a non-mean-field theory, 
taking into account atomic field fluctuations (see Sec. IIIII below). 

In experiments with a trapped BEG the kinetic energy terms in Eqs. pi.20a|) and pl.20bp . arisin g fro m Ha and Hm 
[see Eq. p.3|) ]. can be generally neglected according to the Thomas- Fermi approximation (see Ref. [l02l |). In this case 
(fo (r, t) and (pm (r, t) depend on r only parametrically and Eqs. pi.20a|) and pi.20bp are reduced to a set of ordinary 
differential equations in t. A different situation takes place in the case of an expanding BEG considered below in Sec. 

Eel 

When the kinetic energy terms are neglected, the system can be described by a set of real equations, similar to the 
optical Bloch equations (see Ref. [i^) for the new real variables, the atomic and molecular "populations" (condensate 
densities) 

no(r,i) = |^o(r,^)|^ n„, (r, i) - (r, i) j^, (11.22) 

and "coherencies" 

u (r, t) = 2Re {gifl (r, t) (r, t)) , v (r, t) = -2Im {gipj (r, t) if*^ (r, t)) . (11.23) 

Equations pi.20ap and pl.20b|) lead to the following equations of motion for the new variables, 

no{r,t) = 2v{r,t) -2Ta{r,t)nQ{r,t) (IL24a) 

nm{r,t) = -v{r,t)^2Tm{r,t)nm{r,t) (II.24b) 

v{r,t) ^ D{r,t)u{r,t)~[2Ta{r,t)+Trn{r,t)]v{r,t) 

+2\gfno (r, t) [4n™ (r, t) - no (r, t)] (II.24c) 

u (r, t) = -D (r, t) V (r, t) - [2Ta (r, t) + T^ (r, t)] u (r, t) . (II.24d) 



Here 



and 



D (r, t) = 2e, (t) + V{r) + {2Ua - Uam + 5a) (r, t) 

{r,t) (11.25) 
V (r) = 2Va (r) - Kr, (r) , 



Ta{r,t) ^ ^kanm{r,t) , Tm ir,t) ^ ^kaUo {r,t) + kmfim ir,t) . (11.26) 



A form of such equations similar to the vector form of the optical Bloch equations is derived in Ref. [7£ 
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C. Expanding Bose-Einstein condensate 



Consider an atomic BEC confined in a harmonic trap potential 

3 

2 



2 2 



(11.27) 



In the Thomas-Fermi regime, when the kinetic energy is neghgible compared to the c hemical potential eo, the atomic 
mean field can be described by the stationary Thomas- Fermi solution (see Ref. |102{ ) 



1/2 



ipTF (r) = iVg 



1/2 



^pcak 



TO X ^ 



2 2 



; ^pcak 



15TO^u;ia;2i^3^o ^ 
1287r2x/2^ 



-3/5 



(11.28) 



where iVo is the initial number of atoms, npoak is the peak density, and eo = f^a'^-pcak- 

An expanding BEC is formed when the ato mic tr ap containing the BEC is turned off. The following expansion of 
pure atomic BEC has been considered in Refs. |l03l |. In this case a solution of the single Gross-Pitaevskii equation 



iipo (r,t) 



^0 ir,t) 



can be represented in the form 



using the scaled coordinates 



a scaling factor 



and the phase 



(^0 {r,t) = A{t) $o(p,t) e'^. 



^(r,i)=m^r|M]_ + 5„(t) 



2hj (t) 



(11.29) 



(11.30) 



(11.31) 



(11.32) 



(11.33) 



which incorporates most of the contribution of the kinetic energy. 

The initial conditions at the start of the expansion t = toxp are hj (tcxp) = 1, bj (tcxp) = 0, 5*0 (toxp) = 0, and 
$0 {p,tcxp) — fo (p, ^cxp) — fTF [p)- Substitution of Eq. pi.30p into Eq. (jll.29[) leads to the following equation for 
the transformed mean field $o (p, t) 



-y- 



1 92 



2to ^ 62 Q 2 

1—1 J 'J 



Sq 



*o(p,t). 



(11.34) 



Following Ref. |103j let us take the scales hj as solutions of the set of equations 

bj{t)^uj]A\t)/hj{t) (11.35) 

and Sq (t) = —eoA^ (t). As has been shown in Ref. [l03| | the residual kinetic energy terms can be neglected in the 
Thomas- Fermi regime and Eq. pi.34D is satisfied by the stationary Thomas- Fermi solution $0 {p, t) — ipTF (p)- 

Solutions of Eq. (jll.35|) (see Ref. [lOSj ) demonstrate that the expansion is ballistic after acceleration during a time 
interval t ^ min (w"^). In experiments the expansion is started at a large detuning, when the molecular occupation is 
negligibly small. The resonance is approached and the molecules are formed when the density is substantially reduced. 
Therefore the terms proportional to Um, Uam, Sa, 6m, and Vm (r) in Eq. pi.20p can be neglected. Since the molecules 
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are formed at the ballistic stage of the expansion and inherit the velocity of the atoms they are formed from, the 
molecular field can be represented in the form 



ipm (r, t)^A (t) (p, t) e 



2iS 



(11.36) 



Substitution of Eqs. pi.30[) and pi.36p into Eq. PL20D leads to the following set of coupled equations for the trans- 
formed mean fields (see Ref. 531): 



{t)--A^ [t) K\<^^ {p,t)\^ 



^o{p,t) + 2A{t)g*% ip,t) (p,i) 



-A' (t) f/Jc&o (P, P - 60 + y E ^1^' (P, t) 



AA^ (t) 



fcal^O (P, <)P+fcm|$m (p, t) \ 



$™ {p,t) + A{t) g$^(p,t) 



(11.37) 



(11.38) 



The residual kinetic energy terms can be neglected in the Thomas-Fermi regime as well as in the case of a pure atomic 
BEC. The loss processes and atom-molecule transitions distort $o (p, t) from the Thomas-Fermi shapes, leading to 
additional energy shift compared to the pure atomic case [the last term in Eq. pi.37|) ]. This shift, as well as the 
one expressed by the last term in Eq. pi.38p . is, however, of the order of A^ (t) eo, and can only lead to a negligibly 
small shift of the resonance. Thus, an analysis of an expanding hybrid atom-molecule BEC can be reduced to the 
solution of a set of ordinary differential equations pi.37|) and pi.38p . The parametric dependence on p arises from 
the inhomogeneous Thomas- Fermi initial conditions for $o (p, t) expressed by Eq. pi.28p . 



III. PARAMETRIC APPROXIMATION 



A. Dump state elimination 

Consider an atom- molecule homogeneous quantum gas described by the second-quantized Hamiltonian p.l7p . Let 
the initial state of the atomic field at t = to he a. coherent state of zero kinetic energy 

4>a{p,to)H^{27rf\o{to)6{p)\m), (III.l) 

where \(po (to) = Ua (to) is the initial atomic density and |in) is the time- independent state vector in the Heisenberg 
picture. A pair of condensate atoms forms a molecule of zero kinetic energy. Therefore the resonant molecules can 
be represented by a mean field ipm (t) as 

(inl^-™ (p, t) |in) = {27: ip^ (t) S (p) , (IIL2) 

where \(prn (i) P = (t) is the molecular condensate density. This approach therefore takes into account the time 
dependence of the molecular mean field, but neglects fluctuations of the molecular field due to a Feshbach coupling 
of non-condensate atoms. 

The outcome of atom-molecule and molecule-molecule deactivating collisions is introduced, as in the mean field 
theory of the previous section, by adding the molecular dump states. The elimination of these states in a second- 
quantized description should, however, be done in a diffe rent way (see Ref. [l^). It is similar to the Heisenberg- 
Langevin formalism of quantum optics (see Refs. [I04lll05j ). but takes into account the nonlinearity of the coUisional 
damping. 

The Hamiltonian p.l7p yields the following equations of motion for the the annihilation operators 'i'a (p, t) of the 
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atomic field and 'i'd (p, t) of tlie molecular dump states, 



fa (P,t) 



+2 (27r) 



-3/2 



d^dh (P - p') #1 (p', t) (p + p', t) 



+ (27r) J d^PdPpdd*d (Pd - 2p') ¥^ (p' + Pd - p, t) (p', t) (pd, t) 

d 

(27r)"^/^ ^ y" d^pdd^pmdd (Pd - 2p) *J (pd, t) *m (p,„, i) *a (P + Pd - Pm, t) 



(IIL3) 



i^d{pd,t) = 



Am 



Ed 



"^d (Prf,t) 



(27r)- 



-3/2 



d^p'd^Pmdd (Pd - 2p') (p', t) (p™, i) *a (p' + Pd - Pm, t) 



+ (27r) ^ / d^Pmd^Pudud (P« - Pd) ^'1 (P«> ^) *m (Pm, (P« + Pd - Pm, 

(IIL4) 

The atom and the molecule emerging from the deactivation event p.7|l depart with momenta pd > Pd [see Eq. 
(|IL12P ] . The deactivation energy iJ^j substantially exceeds characteristic energies of atoms formed by dissociation of 
the condensate molecules, allowing us to discriminate two groups of atoms, with momenta above and below min (Pd), 
respectively. [This assumption is equivalent to the condition (jll.3p .] Equations pil.3[) and pil.4[) give the following 
equation of motion for the product of the field operators (with pd > min (Pd)) 



d 



^d (Pd,0*a (P-Pd,i) 



Pd , (Prf-P) 



Ed 



4m 2m 
dd (3pd - 2p) (p,<) iPm (t) 



^-d (Pd,t) *a (p-Pd,t) 



(IIL5) 



where €a is neglected as small compared to Ed, and the molecular field operator (p, t) is replaced by the mean field 
<Pm (t) [see Eq. pll.2|) ]. The source term in Eq. pII.SP arises from the commutation of field operators upon normal 
ordering, while the terms containing the dump field operators are neglected here. Substitution of the solution of Eq. 
piI.SP and the molecular mean field pil.2p into Eq. pil.3[) gives the following integro-differential equation 



i*a(p,i) 



2m 



^a(p,i) + 24 (2p) ip„,it) *t (_p,i) 



lifl, (t) / dt'K (t - t') (^,„ {t') i>a (P, t') + iF (p, t) 



(III.6) 



with a kernel 



Pd_ ^ iVd - p) 
4m 



2m 



K{t-t')=J2j d^'pdldd (3pd - 2p) |2 exp 

d 

and a quantum noise source 

F (p, t) = -lifl, [t) / d^pddd (3pd - 2p) *d (Pd, to) *a (p - Pd, to 



Ed\{t- t') 



X exp 



Pd , (Pd-P) 



4m 



2m 



Ed]it- to) 



(111.7) 



{IIU 



As in the Heisenberg-Langevin formalism, commutators of the quantum noise are related to the kernel by a fluctuation- 
dissipation theorem, except that here this relation involves averages of the commutators. 



F (p, t) , Pt (p', t') jin) = ip*^ (t) ^„ (f) K{t- t') 5 (p - p') 



(III.9) 
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In the Markovian approximation, the kernel is assumed to be sharply peaked a,t t = t', so that 



(III.IO) 



where the deactivation rate coefficient ka is defined by Eq. (jll.lSp . An expression of the kernel by Eq. piI.lOp implies 
that the system retains no memory of its history. The equation of motion for the atomic field then attains the 
Heisenberg-Langevin form 



(p,t) 



Zm 2 



(p, t) + 2g*^,n (t) #1 (-P, t) + iF (p, t) 



(III.U) 



where dh {p) is replaced by its maximal value d/,. (0) = g [see Eq. (jll.6|) ]. The shift associated with the deactivating 
collisions Sa [see Eq. (jll.lGp ] can be neglected compared to other energy scales in real physical situations. The quantum 
noise source F {p,t) is (S-correlated in the Markovian approximation, obeying 



(in| \f (p, t) , Ft (p', t')] |in) = (t) \H {t - t') 5 (p - p') . 



(III.12) 



The Markovian approximation is applied only to deactivating collisions, while the description of the association- 
dissociation processes remains non-Markovian. 



B. Atomic field representation 

Equation pll.lip is a linear inhomogeneous operator equation. Consider at first a solution of the corresponding 
homogeneous equation. Linear operator equations can be solved as c-number ones. A general solution can be expressed 
as 



*a (P, t)=C {t) * (p, to) {P, t) + i'^ (-P, to) {P, t) 



(III.13) 



in terms of the creation and annihilation operators at the initial time t — to and the solutions ipc,s {p, t) of the 
corresponding c-number equations 



2m 



i^c,s {P, t) + 2g*ip^ {t) ipl^^ {p, t) , 



(III. 14) 



given the initial conditions il>cip,to) — 1, V's(p, ^o) = 0- Equation piI.14p as well as the functions ipc,sip,t) are 
independent of the p direction. The factor C {t) takes into account the imaginary term in Eq. pil.lip . 



C {t) = exp I - j dti^\ipm (h 



(III.15) 



A general method for solving inhomogeneous equations is the variation of constants in the solution of the corresponding 
homogeneous equations. Following this method, let us replace the field operators at t = to: which play the role of 
constants in the solution pil.isp . by unknown time-dependent operators A (p, t), representing the atomic field operator 
in the form 



*a (P, t)=C (t) [i (p, t) {p, t) + it (_p^ t) {P, t) 

The initial conditions at t = to require that 

i(p,to) = *a(p,to)- 
Substitution of Eq. (|III.16P into Eq. pll.lip leads to the equation 

i (p, t) {p, t) + it (-P, t) (p, t)^F (p, t) (t) , 



(III.16) 
(III.17) 
(III.18) 



which, together with its hermitian conjugate, form a system of linear equations for A(p,t) and (— p,t) with the 
determinate 



1^, (p,t)|2-|V.(p,t)|^ 



1, 



(III.19) 
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as one can prove by using of Eq. piI.14p . As a result, the operator A (p, t) can be expressed in the form 



t 



/T I 
{P, h) F (p, h) - V'. (p, h) Ft (-P, ti)] . (III.20) 

to 

The reactions (|I.7p of deactivating coUisions thus contribute to the description of the atomic field operator both the 
factor C {t) in Eq. piI.16|) and the second term, containing the quantum noise, in Eq. pil.20p . The factor C {t) 
describes the decay due to deactivating collisions, while the quantum noise provides for maintaining the correct 
commutation relations of the atomic field operators (in an average sense) as 

i!a (p, t) , ^1 (p', t)\ |in) = (5 (p - p') . (III.21) 

C. The atomic mean field and correlation functions 

Physical observables are expressed in terms of averages of the field operators and their products. The atomic field 
operator includes a contribution proportional to the quantum noise F (p, i) [see Eqs. pil.ll[) and pil.20[) ]. Since 
the initial state (described by the state vector |in)) does not contain the molecules A2 (d) and atoms with momenta 
p > Pd, the definition pll.8|l of F (p, t) leads to the following zero-valued averages 

(in|F (p, t) |in) = (in|F (p, t) F (p', t') |in) = (in|Ft (p, t) F (p', t') |in) = 0. (III.22) 

Averages of products of the quantum noise and atomic field operators are zero- valued as well. 
The commutation relation piI.12p thus gives 

(in|F (p, t) Ft (p', t') |in) = ka\fm (t) \H {t - t') S (p - p') . (III.23) 

The average of the atomic field operator is determined by Eqs. (IIII.ip . pll.20|l . and (IIII.16P as 

(inl*^ (p, t) |in) = (271 ipo it) S (p) , (III.24) 

where 

(^0 {t) = C (t) (0, t) (^0 (to) + iJs (0, t) ip* (to)] (III.25) 

is the atomic condensate mean field. 

The averages of the quantum noise lead to the following expressions for the two-atom correlation functions 

(inl^-t (p, t) (p', t) |in) = {2n f no (t) S (p) S (p') + (p, t) S {p - p') (III.26) 

(in|*„ (p, t) i'a (p', t) |in) = (2^)' mo (t) S (p) 6 (p') + {p, t) 6 (p + p') , (III.27) 

where 

no (i) = |(^o (t) P (111.28) 

is the condensate density. 

Us {p, t) = {p, t) P [1 + Vs {p, t)] + IV-c (p, t) \% (p, t) - 2Re [C (p, t) i,, {p, t) {p, t)] (III.29) 
is the momentum spectrum of the non-condensate atoms, and 

IT^oit) = flit) 



(III.30) 



nis {p, t) = 1/;, (p, t) (p, t) [1 + 2rj, {p, t)] - (p, t) 77, (p, t) - (p, t) 7^1 {p, t) 
are the anomalous densities of the condensate and non-condensate atoms. The functions 



t 

to 



(III.31) 



t 

P fit' 

vAp,t) = KC' it) J ^T^l^™ (i')PV's(p,0^c*b,0 

to 
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describe the contribution of quantum noise. 
The atomic density 

na{t) = (27r)"^y dVrf^P2exp[i(p2-pi)r](in|*l(pi,t)#a(P2,t)|in) (III.32) 

then appears to be r-independent, and comprises the sum 

Ha (t) no (t) + n, (t) (III.33) 

of the densities of condensate atoms uq (t) [see Eq. piL28l) ]. and of non-condensate (entangled) atoms Ug (t) in a wide 
spectrum of kinetic energies E — p'^ / (2m), 

ns{t) = J dEfisiE,t), (III.34) 

where the energy spectrum hg (E, t) is related to the momentum spectrum Ug {p, t) [see Eq. (jIII.29p ] as 

fis {E, t) = ^rig (p, t) . (III.35) 

D. The molecular field 

The equation of motion for the molecular field operator (p, t) is obtained by a procedure of dump field elimination 
which is similar to the one presented in Sec. IIII Al for the atomic field. After substitution of the solution of Eq. (|III.5p . 
and a similar equation for the product {Pd,t) ^'m (p„,t), the operator equation of motion attains the form: 

2 r 
i'^m (Pm, t) = (Pm, t) + (2^)"^/' / <f pdh (2p - p^) #a (p, t) 4-, (p^ - p, t) 

t 

~i (2^)-^/' y" dt' J d'^pK {t - t') i!\ (p - p„„ i) if a (P, t') it') 



t 

i {2nf'^ [ dt'Km (t - t') (t) ifl (<') S (p„0 + iFa (p, t) + iF,n (p, t) , (III.36) 



to 



where the kernel Km {t — t') pertains to molecule-molecule deactivating collisions (jl.Sp . In the Markovian approxima- 
tion it is approximated by Km {t — t') = kmS {t ~ t') where the deactivation rate km is defined by Eq. (jII.18[) . The 
quantum noise sources Fa (p,t) and Fm (p,^), related to atom- molecule and molecule- molecule deactivation, vanish 
on mean-field averaging like F (p,t) [see Eq. pil.22p ]. The right-hand side of the resulting equation is proportional 
to 6 (Pm) by virtue of Eqs. (|III.26P and (jIII.27p . thus securing the consistency of Eq. (jIII.2p . Taking into account 
Eq. (jlll.lOp we obtain the equation of motion for the molecular mean field 

oo 

i'pTn {t) = gmo (t) + ^ j p^dpdh (2p) {p, t) - i '^^^ ' (III.37) 



where the the anomalous densities of the condensate and non-condensate atoms, toq (t) and {p,t), are defined 
by Eq. pil.30p and the atomic density is defined by Eq. pil.33[) . These equations can describe the spontaneous 
dissociation of molecules into the atomic vacuum (mg — 0) due to the term in Eq. pil.37[) proportional to iris . 

Equation (IIII.37|) cannot be solved with a constant dh (p) due to the divergence of the integral over p. It can be 
seen from asymptotic behavior of the anomalous density to^ (p, t) as p — > oo. In this limit the solutions of Eq. piI.14p 
can be approximated as 

ipc {p, t) ^ exp ( {t - ^o) 

(III.38) 

2mg* " 



i^s (p, t) 



n2 



ipm (t) exp I {t -to)] - fm (to) exp ( [t - to) 

' 2m / \ 2m 
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leading to 



rris (p, t) 



2mg* 
-— -2— </5m [t) 



and hence to an asymptotically non-decreasing integrand in Eq. (|III.37|) . Following Ref. [i^ let us introduce a 
momentum cut-off Pc, setting dh (2p) = g for p < pc and dh (2p) = otherwise, where g is defined by Eq. (jll.2ip . The 
integral term in Eq. (|III.37p tends to Em^m m the limit Pc oo with — —■^\g\'^Pc- Substituting renormalized 
functions 

{t) ^ (fm (t) CXp {-iE„it) , Vc,s (P, t) ^ Vc.s (p, t) exp (^^Emt 

into Eqs. (|III.14|1 and (|III.37|1 one obtains a non-divergent equation for the molecular mean field 



iifim (t) = giTlQ (t) 



27r2 



dp [p^gm,s {p, t) + 2TO|g| {t)] 



(III.39) 



while Eq. (jIII.14[) retains its form. A numerical solution of Eqs. (jIII.14[) on a grid of values of p, combined with Eq. 
(lfIL39p , is consistently sufhcient for elucidating the dynamics of the system. 

Numerical calculations and a qualitative analysis can be made more convenient by using dimensionless variables. 
Let the time and momentum be rescaled to 



(III.40) 



where the density scale 



n = Ua {to) + 2n„, (io) 



(III.41) 



is determined in terms of the initial atomic and molecular densities. The scale for all energies is thus |(7|\Ai- The 
deactivation rate coefficients are rescaled to 



m 

The equation of motion for the rescaled molecular mean field attains the form 



(III.42) 



i^m (i) = mo (i) + cr / dp [p'^rhs {p, t) + 2(^„ [t]] 



^flQ (t) + I dpp^fls {p,i) + fcm|(pm {t) 1^ 



<^m it) , 



(III.43) 



where the rescaled densities no, fis, i^o^ and rhs are determined by the same equations as the non-rescaled ones, but 
expressed in terms of rescaled parameters. The dimensionless coefficient a is given by 



{m\g\) 



3/2 



27r2ni/4 



{\aa\m^lA)' 



-,1/4 



(III.44) 



expressed in terms of the resonance strength and initial densities. The contribution of non-condensate atoms (atomic 
field fluctuations) to Eq. (|III.43p is proportional to a and, therefore, increases with the resonance strength and is 
suppressed by high densities. In contrast, the role of fluctuations in a pure atomic gas increases with the density, 
according to the Bogoliubov theory (see Ref. 
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E. Relation to the Hartree-Fock-Bogoliubov approach 



The parametric approximation takes into account second order correlations of the atomic field. The Hartree-Fock- 
Bogoliubov approach of Ref. [i^ takes into account the atomic field correlations to the same order, but neglects 
the relaxation due to deactivating collisions. It consists of a solution of the equations of motion for the atomic 
and molecular condensate mean fields (pa it) and Lpm (t), as well as the normal and anomalous densities of the non- 
condensate atoms, Us (j), t) and rris (p, t). A conventional derivation of the Hartree-Fock-Bogoliubov equations is based 
on the Heisenberg equation of motion for the atomic field operator. Use of the Heisenberg-Langevin equation (jllLlip 
allows to include the relaxation process. The definitions of (pa{t), ns{p,t), and ms{p,t) [Eqs. (jIII.25p . pil.29[) . 
and pil.30p . respectively] in a combination with Eqs. (|III.14p and (jIII.3ip lead, after some algebra, to the following 
equations of motion, 



i(po {t) = 



Ea {t) - -ka\(p„i (t) P 



ipa (t) + 2g*(pl (t) <y9„ {t) 



ris {p, t) = -4Im [gip*^ (t) (p, t)] - ka\<f„r (t) {p, t) 



(IH.45) 



ims (p, t) 



rp2 



2ea (t) ~ ika\fm {t) 



rus (p, t) + 2g*(pm (t) [I + 2ns (p, t)] 



When the deactivating collisions are neglected (fc^ = fc^ = 0), these equations, in combination with pil.37p are just 
the momentum representation of the Hartree-Fock-Bogoliubov equations used in [43.] and nothing else, leaving the 
present approach mathematically equivalent to the Hartree-Fock-Bogoliubov one. In the absence of relaxation we 
have rjcs [p, t) ^ and the normal and anomalous densities are expressed as 



m-s [p, t) = {P, t) ipc [p, t) . 



(IH.46) 



Therefore a solution of equations for ipc^^ {p, t) has some benefits for numerical calculations due to less variation of 
these functions compared to Ug {p,t) and {p,t). Moreover, Eq. piI.14p is much easier to analyze, and can even 
produce analytical solutions in some situations. 



F. Exactly soluble models 

Consider the dissociation of a molecular BEC neglecting deactivating collisions {ka = fc^ = 0) and the depletion of 
the molecular field in comparison to its initial value, so that one can approximately assume tpm =const. In this case 
the various atomic p-modes become decoupled and Eq. pil.lip can be solved separately for each mode. The resulting 
equations 

i^a (P, t) = ep (t) (p, t) + 2.r*l (-P, t) , (IH.47) 

2 

where (t) = -\rta {t) and g — gipl^: are similar to the equations for non-adiabatic transitions in two-state systems 
(see Refs. [lOi [iM [iQi)- However, here we have a transition between the atom creation operators ^'jj (— p, t) 
and the annihilation (or hole-creation) one (p, 0- Thus, a transition from a hole to an atom corresponds to the 
formation of an atomic pair as a result of molecular dissociation. Since the field operators are non-hermitian, the 
many-body non-adiabatic transition problem (in contrast to a single-body problem involving a single molecule) has 
an antihermitian coupling matrix, leading to complex adiabatic energies. 

The atomic field operator can be expressed in the form piI.13p with C — 1. The equations of motion for the 
coefficients ipc,s {p, t) can be represented as 

zV-c.a (P, t) = ep it) V'c. (p, t) + 25*V:,c {p. t) , (111.48) 
leaving the initial conditions unchanged, 

V'c(p,io) = l, ^s{p,tn) = 0. (III.49) 

The first-order coupled equations pil.48p for V'c (p, t) and (ps {p, t) can be transformed to second-order uncoupled 
equations 

{p, t) + {el {t) - 4l5l2 + (t)) ^,,s {p, t) = 0, (IH.50) 
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with the same initial conditions (|III.49|) for the two solutions, augmented by initial conditions for the derivatives, 

(to) = -«ep (to) , (to) = -2i~g*, (111.51) 

derivable from Eqs. pil.48p at i = to- The solutions ijjc {p,t) and ips {p,t) can be expressed in terms of a pair of 
conventionally selected independent solutions of Eq. piI.SOp tpi (t) and ip2 (t) , treated as "standard solutions" , in the 
form 

V'c W = WWM [(' ^2 (to) - ti^2 (to)) V-i (t) - (e (io) V-i (to) - iV'i (to)) ^-2 (t) 

^« = w//^*/ 1 (to) V'l (t) - ^1 (to) ^-2 (t)] , (III.52) 

where 

W{V'l,V2}=V'l^2-V'2^1 (111.53) 

is the Wronskian of the standard solutions. 

Consider at first a case of a time-independent energy ep =const. The standard solutions are ipi^2 (t) = exp {±Kt) 
with 



/4|ffP - el, (III.54) 

and ipc.s attain the form (see Ref. [75 

V'c {p, t) — coshK {t — to) — i— sinh n{t — to) , V's {P: t) — —i— sinh K{t — to) ■ (III. 55) 



The modes with ep < 2|g|, or real k, are unstable and grow exponentially at sufficiently long times. This growth 
is related to Bose enhanced dissociation and is similar to the amplification of spontaneous radiation in laser-active 
media, where the parameter k plays the role of the amplification coefficient. A similar solution has been obtained for 
a parametric oscillator in quantum optics, see Ref. 105]. 

If a time-dependent Cp (t) has a real zero, the energies of the hole and atom cross and we face a many-body curve- 
crossing problem (an analog of a single-body curve crossing involving a single molecule, see Refs. [l06lll07Lll08l . ll09l |). 
The simplest curve crossing model is the linear one with a time dependence ep (t) = —f3t, /3 > (see Ref. [83]). The 
second-order equation pil.50p attains in this case the form of a parabolic cylinder equation (see Ref. [llOj ]) 



V^e.s {p, t) + {pH'' - 4l5l2 - zf3) V-,^, (p, t) = 0, (III.56) 

and the standard solutions are expressed in terms of the parabolic cylinder function (see Ref. [87|]). In the limit 
to ^ — 00, t ^ 00, the number of non-condensate atoms formed by the spontaneous process can be written as 

n,(p,t)« (l-e-2-A)e2.A^ 

where 

The parameter A is nothing else but the Landau-Zener exponent (see Refs. [l06l Il07l llOSl [l09| ). since the coupling 
of states, according to Eq. pil.47p . is 2g and the slope of the two-atom energy is 2(3. 

The first factor in Eq. pil.57p is the familiar Landau-Zener probability, while the second one describes an amplifica- 
tion of the spontaneous dissociation due to Bose-enhancement. As in Eq. pil.55[) this amplification is reminiscent of 
lasing, and the exponent can be interpreted as a product of a characteristic crossing time {g / P) and an "amplification 
coefficient" g. 

A distinction between Eq. (IIII.57|) and the Landau-Zener probability can be clarified by interpretation of the crossing 
as scattering on a parabolic potential barrier. Indeed Eq. pII.56P can be represented in the form of one-dimensional 
stationary Schrodinger equation 

- V^e,s (P, t) - iP'^t' - V^e.s {P, t) = -41512^-,,, (p, t) , (III.58) 

where the time plays a role of a coordinate and the "energy" — 4]^]^ is negative and lies below the the barrier (see 
Fig. [2^) . The sign of the energy is a consequence of antihermiticity of the coupling matrix [see discussion following to 
Eq. pil.47p ]. A single-body curve crossing problem has a hermitian coupling matrix. It can be described by the same 
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FIG. 2: Scattering on a parabolic potential barrier corresponding to the many-body crossing (a, barrier penetration) vs. 
single-body crossing (b, above-barrier reflection). The dashed line marks the "collision" energy. 

Schrodinger equation piI.SSp but the energy has an opposite sign and hes above the barrier (see Fig. [^b) . Therefore 
the many-body problem corresponds to barrier penetration with exponentially small penetration probability and large 
reflection in contrast to large penetration and exponentially small reflection corresponding to single-body crossing. 
This leads to the replacement of the negative exponents by positive ones, and yields Eq. pfl.57p in place of the 
Landau-Zener formula. 

The method of Ref. [13] has been applied to a case of heteronuclear molecules in Ref. (s^l. Exact sol ution s for 
some other shapes of Cq (t) have been obtained using the theory of quantum integrals of motion (see Refs. jlll| and 
references therein). 

G. Bose-enhancement and normal density 

The effect of Bose enhancement is mathematically expressed as exponentially growing solutions of the equations 
of motion for the atomic field. Due to equivalence of the parametric approximation and Hartree-Fock-Bogoliubov 
approach, Eqs. (|III.45P have in this case exponentially growing solutions as well. With the neglect of deactivating 
collisions and of depletion of the molecular field, and for £p =const, the equations for the anomalous and normal 
densities attain the form of a system of inhomogeneous linear differential equations, 

im, {p, t) = 2epm, (p, t) + 2g* [1 + 2n, [p, t)] (11159) 
(p, t) = 2~g*ml (p, t) ~ 2gm, (p, t) . (III.60) 

The corresponding homogeneous system has solutions exp (±Kt). The unstable modes with < 2|g|, or real n [see 
Eq. (|III.54p ] demonstrate exponential growth. 

Some approximate approaches are based on the neglect of the normal density (p, t) — Q. In this case Eq. (|III.59p 
has only an oscillating solution ~ exp (±2iepi) and can not describe the effect of Bose enhancement. As a result 
such approximations met little success in some situations (see Ref. Ill2j |). The normal density have to be taken 
into account also for a correct description of the squeezing of non-condensate atoms (see Se c. IVIII . Nevertheless, 
the justified neglect of the normal density in the first-order cummulant approach of Refs. [13, l83l. |84|. [ssl . [s^ leads 
to significant simplification of resulting equations and allows for the description of effects of spatial inhomogeneity 
beyond the local density approximation. The normal density could be taken into account in higher-order cummulant 
approaches. 

IV. LOSS OF ULTRACOLD ATOMS 

The effect of Feshbach resonance in ultracold atomic collisions has been observed at first accompanied by a dramatic 
loss of trap population when the magnetic field was tuned to the vicini ty of the resonance. Such loss has been observed 
both with atomic BEC [H, HI, [s^ and with ultracold thermal gases (ll3j | . 

These experiments can be generally divided to two categories, called in the pioneer MIT works [13, [l^ as " slow 
sweep" and "fast sweep" ones. In the fast sweep experiments the magnetic field has been ramped through the 
resonance, while in the slow sweep ones the ramp has been stopped short of crossing. The experiments of the second 
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kind use generally slower ramp speed than the fast sweep ones, or even a static magnetic field, although the principal 
difference between the two kinds of experiments is the presence or absence of resonance crossing. 

The atomic loss is mainly attributed to two loss mechanisms. The first one is a collision-induced deactivation 
process, such as (II. 7[) and p.Sp (2^ . [stI [40j . undergone by the temporarily formed molecular BEC. The second 
mechanism involves the molecular dissociation induced by crossing to non-condensate atomic states [s^ IH, . The 
effects of the loss mechanisms are non- additive (see Ref. An additional loss can be produced by secondary 

collisions of the relatively hot products of the reactions p.7[) and p.8[) with condensate atoms (see Refs. |37l |45| ) . 



A. Deactivation loss mechanism 



This loss mechanism is included in the mean field approach discussed in Sec. |TT1 Certain properties of the loss 
process can be described by analytical expressions, obtained from Eq. pi.24[) . whenever the following "fast decaif 
conditions hold (see Ref. [40|): 

/iS « r,„ {D + Tl^/D) , r„ » r, (iv. la) 

i?2+rf„»6|5|2no. (IV.lb) 

For definitions consult Eqs. pi.25p and pi.26[) . These conditions mean that the relaxation of n^, v, and u is much 
faster compared to that of no and to the rate of change of the energy, caused by the magnetic field with a ramp speed 
B. Therefore the values of the fast variables can be related to a given value of the atomic condensate density np, 
using a quasi-stationary approximation, by 

and the condition pV.lbjl leads to <JC no. As a result, a single non-linear rate equation for the atomic density can 
be extracted. When terms proportional to the atomic and molecular densities in D [see Eq. pi.25P ] are neglected, the 
resulting rate equation is 

. X Slgpfcan^ (r,i) 

no r, t) = ° ' ^ . IV.3 

[2e,(t) + y(r)]V[fcano(r,<)/2]2 

(The neglected terms in D effectively add an extra shift to the resonance, but its contribution is hardly noticed in 
the present problem.) 

Equation pV.SP has a form anal ogou s to the Breit-Wigner expression for resonant scattering in the limit of zero- 
momentum collisions (see Refs. [39l . lll4| ). In the Breit-Wigner sense one can interpret kafiQ as the width of the decay 
channel, while the width of the input channel is proportional to \g\'^ . This observation establishes a link between the 
macroscopic approach used here and microscopic approaches that treat the loss rate as a collision process. However, 
the right-hand side of Eq. pV.3P is four times smaller than the usual Breit-Wigner expression. A factor of | is 
associated with the loss of a third condensate at om in the reaction p.7[) . while an additional factor of ^ is associated 
with the effects of quantum statistics (see Refs. [I15l Ill6l |). 



1. Approaching the resonance 

Very close to resonance the behavior of Eq. pV.Sp effectively attains a 1-body form linear in n. But as long as we 
stay out of this narrow region, by obeying the "off-resonance^^ condition , 

kan{v,t)^\V{v)~ ^lB{t)l (IV.4) 

we can write Eq. (to a very good approximation) in the 3-body form 

^ +^ 3 27r|aa|fcaA lr\T K\ 

no = -SiCs r, t) ng, = —— — _— _. IV.5 

m/i [B (t) - Bq-V (r) / fi\ 

The three-body rate coefficient is defined here according to conventional chemical notation. It is three times less 
then the one in Refs. [H, [13, SO] • The dependence of Eq. pV.Sp on the scattering length Oa follows from Eq. pi.2ip . 
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FIG. 3: The surviving mean density (see Eq. (|IV.7|I vs. the stopping value of the magnetic field. The resonance was approached 
fi'om below with two ramp speeds, 0.13 G/s and 0.31 G/s, or from above, with the ramp speed -0.06 G/s. These are compared 
with the experimental results fssj] (squares, triangles and circles). 



The above results are obtained neglecting the atomic transport described by the kinetic energy terms in Eq. (jll.20|l . 
This approximation is valid while the characteristic times are small compared to the trap periods. A different situation 
takes place in the slow-sweep experiments 0, (35| , where the characteristic sweep times are large compared to the 
trap periods, and the atomic transport has time to restore the Thomas-Fermi distribution. In this case the averaging 
of Eq. PV.SP over the distribution pi.28|) leads to the following rate equation 



50 

n = -3— i^3(r, t) 



(IV.6) 



for the mean BEC density n (t) — |npoak (t)- 

When the magnetic field ramp is assumed to vary linearly in time, starting at to and ending at t, and Eq. (|IV.4|) 
applies throughout the ramp motion (i.e., by avoiding passage through the resonance as in the slow sweep experiments), 
this rate equation can be solved analytically as, 



n{t) = n{to) 



200 , 9 , , , 



-1-1/2 



(IV.7) 



where B is the magnetic-field ramp speed and the extrapolated time of reaching exact resonance is chosen to be t = 
so that both t and to have the same sign. A similar expression for the homogeneous case has been presented in Refs 



[37|, |40'|. Figure [3] demonstrates a reasonably good agreement between the model and experimental data of Ref. [35 1 
for the 907 G resonance in ^'^Na. The model uses the parameter values A = 0.98 G, Ua = 3.4 nm, and fj, = 3.65/zb, 
calculated in Ref. [3^, and ka — 5.5 x lO^^^cm'^/s, measured in Ref. 12], and does not contain adjustable parameters. 



2. Passing through the resonance 

The off-resonance approximation of Eqs. (|IV.4[) and PV.SP does not hold very close to resonance, and is therefore 
inapplicable to the description of the fast-sweep experiment, in which the Zeeman shift is swept rapidly through the 
resonance, causing dramatic losses. Nevertheless, the fast decay approximation pV.l[) may still be valid. A simple 
analytical expression can then be derived for the condensate loss on passage though the resonance if, in addition, the 
magnetic field variation lasts long enough to reach the asymptotid^ condition 



fi6B > kaUo, 



(IV.8) 
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where 6B is the total change in B accumulated over the sweep. This condition allows the extension of the ramp 
starting and stopping times to =Foo- 

One can now evaluate the variation of uq (r,f) in the infinite time interval (— cx),cx)). Let us rewrite Eq. pV.3p in 
the form 

dnoir^ 3fcano(r,t)|gP 



"0 (^^lBty + ikano{r,t)/2) 



2 



where V (r) is removed by our choice of the origin on the time scale. We then integrate the left-hand side with respect 
to uq from uq (r, — oo) to no (r, oo) and the right-hand side with respect to t from — oo to oo, considering uq as a 
well-defined function of t. The latter integral may be evaluated by using the residue theorem, closing the integration 
contour by an arc of infinite radius in the upper half-plane. The integral along this arc vanishes according to the 
asymptotic behavior of uq (see Ref. '40']). 

The final result does not depend on ka and has the form (valid for all positions r) 

(r. '^o (j^^-oo) 67r|ffp 12TT^\aa\ A mrm^ 
no (r,oo = ; -, s — — — —. (IV. 10) 

^ ' 1 + sno (r, -c3o) ' ^\B\ m \B\ ^ ' 

The product sno in Eq. pV.lOp would be proportional to the Landau-Zener exponent for the transition between the 
condensate and the resonant molecular states whose energies cross due to the time variation of the magnetic field, if 
one could keep the coupling strength gipQ constant. However, for the non-linear curve crossing problem represented 
by Eq. pi.20p . the Landau-Zener formula is replaced by Eq. pV.lOp . which predicts a lower crossing probability, 
since the coupling strength gipo decreases, along with the decrease of the condensate densi ty during the process. 



The transition probability in the Landau-Zener model with relaxation, obtained in Ref. 117|, is given by the same 
Landau-Zener formula, as in the absence of the relaxation, and is independent of the relaxation rate as well as Eq. 
pV.lOp . However, unlike the Landau-Zener formula, Eq. (llV.lOP is applicable only at strong relaxation, when the fast 
decay approximation is valid. 

The asymptotic result pV.lOp describes the decay of the condensate density. Assuming a homogeneous initial 
density within the trap, Eq. pV.10[) applies also to the loss of the total population No{t) = J no(r,t)d^r. An 
asymptotic expression for the total population can also be found when the homogeneous distribution is replaced by 
the Thomas- Fermi one [see Eq. pi.28p ]. In this case, given npcak is the maximum initial density in the center of the 
trap, one obtains 



7Vo(oo) _ 15 / 1 1 1 / 1 



No i-oo) 2snpcak I 3 



\ — . (IV.U) 
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B. Dissociation loss 



The consideration above did not include the decay of the resonant molecular state by transferring atoms to excited 
(discrete) trap states or to higher-lying non-trapped (continuum) states. In a forward sweep, when the molecular state 
crosses the atomic ones upwards (see Fig.^l), this process can be described by curve-crossing theory as a sequence of 
two-state crossings. 

Consider the system in a normalization box of volume V. The number of states of a pair of atoms with momenta 
p and — p per unit of the energy e = /m of their relative motion is 

Due to the time variation of the magnetic field, the molecular state crosses at the time t the atomic pair state with 
the energy e — — 2ea {t) — /i {B (t) — Bq), altogether crossing u^^B such states per unit time. If the depletion of the 
molecular field during each crossing can be neglected, the crossings can be described by the many-body curve-crossing 
model of Sec. IIIIFI The number of atoms formed by each crossing is determined by Eq. (|III.57p . with 

A = TT;]-"™- (IV.13) 

7n\B\ 

Therefore, the loss rate of the molecular population is 

iV„ = -i(l-e-2-A^e2.Ve/iS. (IV.14) 
If the molecular density is small enough, such that A <C 1, one obtains a loss rate of the molecular density 

n,n « -2rdisnm, Tdis = |aa/i| Aa/ {B [t) - Bq), (IV.15) 

which is proportional to the dissociation width Fdis- The same result can be obtained by using of Landau-Zener 
theory (see Ref. [13]). Therefore one can account for the decay of the resonant molecular state into excited trap 
states by adding a term — iFdis^'m to the right hand side of Eq. (jll.20b[l . This approach is in good agreement with the 
parametric approximation only whenever (a) the variation of the molecular field during each crossing is negligible; and 
(b) A <C 1 and the quantum effect of Bose enhancement, described by the positive-exponential factor in Eq. pV.14[) . 
is negli gibl e (see Sees. IIIIFI and IIII Gp . These conditions are obeyed with the parameters used in the calculations 
of Ref. [40[ for the Na resonances , the results of which are confirmed by using the parametric approximation. For 
example. Fig. [5] compares the two approximate methods for the Na resonance at 853 G with the strength A « 9.5 
mG (see Ref. [39], other parameters having the same values as for the 907 G resonance presented at the end of Sec. 
IIVAIII . 

However, under other conditions the two methods produce different results. Figure [6] compares them for a case of 
^^Rb resonance (see Ref. 56]). The values of A = 11 G, Bq — 154.9 G, and Qa — —450 (in atomic units) are taken from 
Ref. [lis], and the value of /i = — 2.23/iB is taken from Ref. [131 • The condensate loss calculated with the parametric 
approximation is more then two times higher then the one calculated with the mean field approach. The main 
loss mechanism is the Bose-enhanced crossing to the non-condensate atomic states. The Bose-enhanced dissociation 
drastically reduces the molecular condensate density, compared to the mean field results, leading to the suppression 
of deactivation losses (the total density of condensate and non-condensate atoms remaining almost unchanged and 
the results are insensitive to deactivation rates for ka,km < 10~^cm'^/s). As a result, the parametric approximation 
demonstrates a much better agreement with the experimental data than the mean field results (see Fig. [7|). The 
drastic increase of loss at slower sweeps is actually related to the effect of Bose-enhancement, as demonstrated by the 
plot for the maximal value of the non-condensate state occupation Ug {p, t) in Fig. [71 The disagreement between the 
parametric approximation and the experimental data at intermediate ramp speeds can be related to fluctuations of 
the molecular field and higher-order correlations of the atomic field neglected in both theories. 

It is interesting that the inclusion of the dissociation loss mechanism reduces the total loss (cf. the dotted and dot- 
dashed lines in Fig.. [7]). This paradoxical result has the following explanation. The coupling to non-condensate atom 
continuum shifts the molecular state energy (above the continuum threshold the shift is imaginary and proportional 
to the width). The energy-dependent shift effectvely accelerates the sweep, reducing the probability of crossing from 
the atomic to molecular BEG. For weaker resonances jwhen the crossing probability remains almost unchanged, an 
additional effect becomes most prominent (see Ref. 001 )• The atom- molecule deactivation p.7p leads to the loss 
of three condensate atoms per each resonant molecule formed, while the dissociation leads to the loss of only two 
condensate atoms. Therefore, the inclusion of dissociation losses transfers flow from the deactivation (11.71) to the 
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FIG. 5: Time dependence of the densities of the atomic (solid lines) and molecular (dashed lines) condensates calculated for 
the Na resonance at 853 G using the parametric approximation (black lines) and the mean-field approach taking into account 
the dissociation width Fdis (red lines). 
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FIG. 6: Time dependence of the densities of the atomic (solid lines) and molecular (dashed lines) condensates calculated 
for the *^Rb resonance using the parametric approximation (black lines) and the mean-field approach taking into account 
the dissociation width Fdis (red lines). The dot-dashed line plots the density of non- condensate atoms calculated using the 
parametric approximation. The ramp speed is B = 5 G/ms. 



dissociation, thus reducing the number of condensate atoms lost per each resonant molecule formed. By the same 
reason in a case of small dissociation loss, the increase of molecule-molecule deactivation rate coefficient reduces the 
total loss (see calculations for Na in Ref. (40|). 



V. FORMATION OF MOLECULES 



In the experiments 0, H, S Hi [IHj [13] diatomic molecules have been formed by sweeping the Zeeman shift through 
resonance in a backward direction, so that the molecular state crossed the atomic ones downwards (see Fig.[S]). This led 
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FIG. 7: Loss of atomic condensate (solid line) and maximal non-condensate state occupation (dashed line) as a function of 
inverse ramp speed calculated using the parametric approximation in comparison with the experimental data by Cornish et 
al. [3^ (circles). The dot-dashed and dotted lines demonstrate, respectively, the results of mean-field calculations taking into 
account the dissociation width Fdis and without this account for fca = 2 x 10~^ cm? /s. 
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FIG. 8: Schematic illustration of transitions between atomic and molecular condensates and non-condensate atoms on a 
backward sweep. 



to the transfer of population from the lowest atomic (condensate) state to the molecular state, as had been proposed 
in Rcf. 39]. In a backward sweep, unlike the forward one (see Fig. Bt , the transitions to non-condensate atomic states 
are forbidden in semiclassical Landau-Zener-type theories (see Refs. jlQ6l . IIOTI llOSl [i09| ) as the second crossing, from 
the molecular to a non-condensate atomic state, precedes in time the first one, from the atomic to the molecular 
condensate state. Nevertheless, quantum theory allows the non-condensate states to be populated temporarily by 
so-called ^^counterintuitive transitions^^ (see Ref. [54{), which are most notable in strong resonances or at low densities 
(see Ref. [111). The molecules are also lost due to deactivating collisions with atoms or other molecules (|I.7p and 
p.8|) . Therefore, an adequate analysis of molecular formation requires a quantum many-body theory that takes into 
account relaxation, such as the parametric approximation. 

Figure [9] compares various approximations for the case of the 1007 G resonance in ®^Rb. The parameter value s 
A — 0.2 G and Oa = 99 atomic units have been measur ed in Ref. [ll9j |. /i = 2.8/iB has been calculated in Ref. [113| |. 
and ka = 7 X 10~^^ cm^/s has been estimated in Ref. [114| . The molecule-molecule deactivation rate coefficient is 
estimated here as km = 10~^''cm'^/s. 
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FIG. 9: Results of calculations for the 1007 G resonance in *^Rb using the parametric approximation with deactivation (black 
lines) and without it (green lines). The red lines plot the mean-field results. Part (a) presents the fractions of the atomic 
populations surviving in the condensate form (solid lines) and converted to non-condensate atoms (dashed lines). The fractions 
of atomic population converted to molecules are plotted by the solid lines in the part (b), while the dashed lines plot the sum 
of all atomic and twice molecular populations. All plots are calculated for the initial atomic density 10^^ cm~'^ and the ramp 
speed of 1 G/s in a backward sweep. 



The results of calculations using the parametric approximation demonstrate that the temporary non-condensate 
atom population persists during a time comparable to the deactivation time. This population substantially reduces 
the maximal molecular density and shifts the maximum time compared to the mean- field results. This differe nce i s 
less pronounced in the case of the weaker 685 G resonance (see Fig. [To]), with A = 17 mG, = 1.4^b (see Ref. fTTsj, 
other parameters are the same as for the 1007 G resonance). 

Figures [9l and fTOl contain also results of parametric calculations neglecting the deactivation losses (fca — km = 0). 
In this case the molecular density increases monotonically in time due to the association of non-condensate atoms, 
while the atomic density remains almost unchanged. 

The conversion efficiency 2max(?i,„ (t)) /no (to) is determined by a concurrence of three processes: the association 
of the atomic condensate and the two loss processes — the dissociation of the molecular BEG onto non-condensate 
atoms and the deactivation by inelastic collisions. A qualitative behavior of the conversion efficiency can be analyzed 
using the rescaled equation Eq. pil.43p . The first term in the right-hand side is related to association, the second one 
— to dissociation, and the rest — to deactivation. The problem depends on the parameter a [see Eq. (IIII.44I) ]. the 
scaled deactivation rate coefficients fca,m, and a scaled ramp speed, which can be expressed as {sn)~^ in terms of the 
parameter s defined in Eq. pV.lOp . Given the rescaled ramp speed, an increase of the resonance strength A increase 
the parameter a [see Eq. (jin.44[) ] and, therefore, the dissociation term in Eq. (jIII.43p . but decrease the deactivation 
rates proportional to ka.m [see Eq. (jin.42p ]. A decrease of the atomic density has the same effect. 

The part of associated atomic condensate increases with decrease of the ramp speed, saturating at values of sn ~ 1 
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[see Eq. pV.10|l and the following discussion]. The molecular loss is proportional to the time and, therefore, it is 
inverse proportional to the ramp speed. This leads to the existence of an optimal ramp speed B oc nA. 

This qualitative analysis is confirmed by results of numerical calculations (see Figs. [Tl1and ll2p . The maximum in 
the ramp-speed dependence is not reproduced by calculations neglecting the deactivation. The mean-field calculations 
predict the maximum, but overestimate the conversion efficiency and provide a slower optimal ramp speed, especially 
for the stronger resonance (see Fig. Ilip . The conversion efficiency at the optimal ramp speed decreases at high 
densities due to deactivation and at low densities due to non-condensate atom population (see Fig. I12p. reaching a 
maximum at the density 10^^ cm~'^ for the weaker 685 G resonance. In the case of the stronger 1007 G, resonance 
the non-condensate atom population is the dominant loss process in the full range of densities considered here, and 
the conversion efficiency remains almost unchanged, slightly increasing with the density. 

The low density required for a more effective conversion can be achieved by the use of expanding condensates, as 
in experiments Si, 113. In the s^Rb experiments the magnetic field ramp was started after a preliminary 

expansion interval tp > 2 ms, following the switching off of the trap. In this case the analysis can be satisfactorily 
carried out by using the mean-field approach (see Ref. [ssj and Sec. Ill C|) . This is demonstrated by Fig.[T31 comparing 
the atomic and molecular densities calculated with the parametric approximation to those calculated with the mean 
field approach, under the appropriate conditions. The initial atomic density used corresponds to the mean density 
reached at the expansion time of 2 ms. As one can see, already when the magnetic field is just 0.5 G below the 
resonance (in a sweep totaling 2.2 G), the results of the two kinds of calculations coincide. When faster magnetic 
sweeps or higher densities are considered, the results of the two approaches converge even faster. 

The applicability of the mean-field approach at the experimental conditions allows us to analyze the case of an 
expanding condensate using a numerical solution of Eqs. pi.37p and (|II.38|) (see Fig.[T31). The expansion reduces the 
atomic density, leading to a slower loss of atoms and molecules. The results are in agreement with the experimental 
data of Ref. [9] reporting that ~ 7% of the atoms are converted to molecules and ~ 30% remain in the atomic 
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FIG. 11: Conversion efficiency as a function of the ramp speed calculated for the 1007 G (circles) and 685 G (triangles) 
resonances in *^Rb with the initial atomic density of 10^^ cm~'^ using the parametric approximation with deactivation (solid 
lines) and without it (dotted lines). The dashed lines plot the mean-field resufis. 
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FIG. 12: (a) conversion efficiency and (b) the lifetime of the molecular condensate Tm (solid lines) at the optimal ramp 
speed {dB/dt)^^^. (dashed lines), all plotted as functions of the initial atomic density no (to), calculated with the parametric 
approximation, for the resonances at 1007 G (circles) and 685 G (triangles) in *^Rb. 
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FIG. 13: Fraction of atomic population surviving in the atomic condensate (a) and converted to molecules (b), calculated using 
the parametric approximation (black lines) and the mean field approach (red lines), for the 1007 G resonance in a homogeneous 
*^Rb BEC with an initial atomic density 2 x lO'^^cm"^ and a magnetic ramp speed of 1 G/ms. The fraction of non-condensate 
atoms calculated with the parametric approximation is represented by the dashed line in part (b). The green lines represent 
the fraction of surviving and converted population calculated using the mean field approach while taking account of spatial 
inhomogeneity and expansion. 



condensate for ramp speeds of less than 2 G/ms. 



VI. FORMATION OF ENTANGLED ATOMS 



The dissociation of the molecular BEC in a forward sweep (see Fig. |4]) has been considered in Sec. IIVI as one of 
the loss mechanisms. However, the dissociation is not just an undesirable process, since the molecules dissociate 
to entangled pairs of atoms with opposite momenta (see Ref. [55*1), as in the case of an unstable atomic BEC (see 
Ref. 89]). Indeed, substitution of the molecular mean field pil.2[) into the atom-molecule coupling p.21|l leads to 
an interaction proportional to a product of atomic field operators with opposite moment a. Interac tions of similar 
form are known in the theory of parametric down conversion in quantum optics (see Refs. jl04 Il05j |). Likewise, the 
atoms are formed in two-mode squeezed states, which are similar to the state of electromagnetic radiation formed by 
parametric down conversion. 

In order to clarify the nature of this entanglement, let us write out a pseudo-Hamiltonian 

Ha=Y.H{p,t), (VLl) 

P^>0 

that leads to equations of motion for the atomic field (|III.11[) . expressed as a sum of contributions of different 
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momentum modes in a normalization box, 

Ha (P, t) 



%i (p, t) i>a (P, i) + *1 (-P, i) *a (-P, t) 



2dh (2p) (i) (p, t) (-P, t) + iF (p, ¥^ (p, + ^-t (-P, t) 



h.c. 



Although the pseudo-Hamiltonian is not hcrmitian, it allows writing the time evolution operator in the form 



P^>0 



where (using the time-ordering operator T) 



17 (p, t) = T exp I j Ha (p, t') dt' 



(VI.2) 



(VI.3) 



The representation of the operator Q, [t] as a product of single- mode operators Cl (p, t) follows from the commutativity 
of the Ha (p, t) with different values of p. 

Imagine a measurement, represented by a projection operator P(p), which selects atoms with the momentum p 
moving in the positive z-direction, and does not affects atoms moving in the negative z-direction. This measurement 
reduces the state vector Q, (t) |in) to P (p) Cl (t) |in) = (l (p, t) |in). The distribution of atoms moving in the negative 
z-direction (the average number of atoms with the momentum p') after this measurement will be determined by 



{inp (t) *t (p', t) (p', t) n (p, t) |in) cx 



(VI.4) 



representing an entanglement of atoms with opposite momenta. This analysis is similar to the one used in the 
entanglem ent of the signal and the idle in the process of degenerate two-photon down-conversion in quantum optics 
(see Refs. |l04l . llOSf ). The state of atoms is also perfectly number-squeezed, i. e. it has zero variance in the relative 
number of atoms with opposite momenta, as in the case considered in Ref . [57| . 

As demonstrated in [87| the non-condensate atoms produced by molecular dissociation are formed in squeezed 
states, which now turn out to be two-mode squeezed states, as in [8^. As in quantum optics, the squeezing is related 
to the quadrature operators 



^a iP,t)±^a i-p,t) 



The uncertainties of the quadratures can be written out as 

1 



{m\X (pi, t) X (p2, i) |in) = -(5 (pi - P2) { 1 + 2n, {p, t) ± 2Re [m, (p, t) e"'"] ] 



(VI.5) 



(VI.6) 



where the momentum spectra (p, t) and rag {p,t) are defined by Eqs. pil.29p and pil.30p . respectively. 

The uncertainties attain maximal and minimal values at two orthogonal values of the phase angle 9. The amount 
of squeezing can be quantified by the energy-dependent parameter 



, 1, (in|X(pi,t)X(p2,t)|in)„ 

r (e, t) = - In 

4 (in|X(pi,i)X(p2,t)|in), 



^ 1 l + 2n, {p,t) + 2\ms{p, t)\ 
4 \l + 2ns{p,t)-2\ms{p,t)W 



A mean squeezing parameter, weighed by the spectral density of Eq. pil.34[) . 



dEhs iE,t)r{E,t)/ns (t) . 



(VI.7) 



(VI.8) 



is used to describe the time variation of the squeezing. 

In the case of a rather weak resonance a description of the dissociation of a pure molecular BEC can be attempted 
by a simple analytical model based on the approach of Sec. IIVBI (see also Refs. [l^, [H, H^). The molecular density 
is given by the solution of Eq. pV.15p , 



n-m (t) = (td) exp 



-2|a,M|A J dt'^mpL {B [t') - B^) 
td 



(VI.9) 
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t (lis) E (mk) 

)KFIG. 14: (a) Fraction of surviving molecular population calculated by using Eq. (|VI.9|) (solid line) and by the parametric 
approximation for the initial molecular density (td) = 10^° cm~^ (dotted line), 10^^ cm~^ (dot-dashed line), and 10^^ 
cm~^ (dashed line), (b) Energy distribution calculated by using Eq. (|VI.10|| (solid line) and the parametric approximation 
for the initial molecular density 10^^ cm~^ (x) and 10^^ cm~^(+). 



where td is the time of start of the dissociating ramp. A molecule dissociates into two entangled atoms, each with an 
energy E, by a crossing occurring at fi [B (t) ~ Bq) = 2E. The resulting energy distribution of the formed atoms can 
be written out as 



ns (E) 



1 



^^^peak exp 



- (E/Epcak)^^'^ 



(VI.IO) 



where 



2m 



m\B\ 
4|aJA 



2/3 



(VI. 11) 



is the peak energy. Although the model does not produce the correct molecular density dynamics (see Fig. [T4k). it 
gives a surprisely good agreement with the numerical results for the energy distribution (see Fi g. | 14b). The analytical 
model demonstrates also a very good agreement to experimental data, as reported in 

Refs.[lOrtll- It should be 

noted that at high energy Eqs. (|III.29p . (|III.35|1, a nd (|III.38P lead to fis (E) cx E'^^^. Therefore, the distribution has 
a divergent mean energy. The peak energy (IVI.ll I) is well defined for each distribution. It is proportional to the mean 
energy for the distribution (|VI.10[) used in Rcfs. [id IT^ . 

In the case of molecular condensate dissociation the squeezing parameter does not exceed the value of r 0.35 
even for the initial molecular density of 10^^ cm~^. The nature of such low squeezing can be understood by using the 
exactly soluble many-body curve-crossing model of Sec. lIIIF] Using expressions for the uncertainties of the quadrature 
operators (|VI.5|) derived in Ref. [s^l, one can represent the squeezing parameter (|VI.7|) as 



2tt\ 



-l)'/' + e-^ 



I (e2^^ - 1)' 



(VI.12) 



where the Landau-Zener parameter A is defined by Eq. (jIV.13|) . In the case of A ^ 1, when Eqs. (jVI.9p and (|VI.10|) are 
applicable [see discussion following to Eq. (jlV.lsp ], the squeezing parameter is r « ^/2nX <C 1. In the opposite case 
of A ^ 1, Eq. (jVI.12[) relates the squeezing parameter to the occupation of the non-condensate mode Ug {p, t) « e^"^^ 
[see Eq. (jIII.57| as r = ^ I n (4ns), i u agreement with an exact relation Ug — sinh^ r, derivable from the definition of 
squeezed states (see Refs. |104l . llOSf ). This relation between high squeezing and high mode occupation demonstrates 
that a correct calculations of the entangled atomic pair properties require taking into account the normal density. 

Higher squeezing can be obtained on dissociation of molecules formed temporarily from the atomic BEC on a forward 
sweep (see Fig.|4|). In this case (see Fig.[T5|) the molecular density is very low (about 10^° cm~^) and persists a shorter 
time (compared to that in the backward sweep) due to fast dissociation, leading to a negligibly small deactivation 
loss. Figure [T5k demonstrates almost full transformation of the atomic BEC into a gas of entangled atomic pairs in 
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FIG. 15: (a) Time dependence of the densities of the atomic condensate (dotted line), entangled atoms (solid line), and the 
molecular condensate (dot-dashed line, xlO'') calculated for the 685 G resonance in *^Rb with the initial atomic density lO'^^ 
cm~^ and a ramp speed of 20 G/s in a forward sweep. The dashed line shows the mean squeezing parameter f [t) [see Eq. 
(|VI.8|l ]. (b) Energy spectra of the entangled-atom density hs {E,t) (solid line) and the squeezing parameter r{E,t) (dashed 
line) calculated at f « 23 ms. 



two-mode squeezed states with the mean squeezing parameter f « 2.4, corresponding to a noise reduction of about 
20 dB. The energy spectra of the entangled-atom density and the squeezing parameter are presented in Fig.fTSb. The 
density spectra are rather narrow, and the peak energy increases with time. The squeezing parameter reaches the 
value of r {E, t) « 3.4 (corresponding to noise reduction of about 30 dB) for the mode with occupation rig « 240 at 
the energy E ^ 1 nK. Comparable squeezing of maxr(i?, i) « 3.6 have been calculated in Ref. ^55^ for the 853 G 
resonance in Na. These values substantially exceed vtvsxr {E,t) « 1.4 calculated in Rcf. 56] for the much stronger 
resonance in ®^Rb. 



Conclusion 



A correct theoretical treatment of Feshbach resonance in BEC has to take into account both quantum fluctuation 
(formation of non-condensate atoms) and damping due to deactivation of the resonant molecules in inelastic atom- 
molecule and molecule-molecule collisions. The deactivating collisions lead to condensate loss and the limitation of the 
atom-molecule conversion efficiency at high densities and slow sweeps. The quantum fluctuations of the atomic field 
describe Bose-enhanced condensate losses, the limitation of atom-molecule conversion efficiency at low densities, and 
the formation of entangled atomic pairs in two-mode squeezed states. The effects of fluctuations and deactivation on 
the loss are non-additive. The squeezing and Bose-enhancement of non-condensate atom formation can be described 
by second-order correlations, but only if both normal and anomalous densities are taken into account. 

The parametric approximation takes into account both the damping and fluctuations in a proper way. This method 
gives analytical solutions for some model systems and can be used for numerical calculations in more general cases. 
The results of numerical calculations generally agree with the experimental data and provide optimal conditions for 
atom-molecule conversion and quantum properties of entangled atoms. 
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